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Abstract 

We resum to next-to-leading order (NLO) the distribution in the light-cone momentum p + = Ex — \px\ and 
the spectrum in the electron energy E e in the semileptonic decays B — > X u lv 1 where Ex and px are the 
total energy and three-momentum of the final hadron state X u respectively. By expanding our formulae, we 
obtain the coefficients of all the infrared logarithms a 7 g L k at 0(a|) and at 0(a|), with the exception of the 
ajl coefficient. By comparing our 0(a 2 s nf) result for the p+ distribution with a recent Feynman diagram 
computation, we obtain an explicit and non trivial verification of our resummation scheme at the two-loop level 
(rif is the number of massless quark flavors). We also discuss the validity of the so-called Brodsky-Lepage- 
Mackenzie (BLM) scheme in the evaluation of the semileptonic spectra by comparing it with our results, finding 
that it is a reasonable approximation. Finally, we show that the long-distance phenomena, such as Fermi motion, 
are expected to have a smaller effect in the electron spectrum than in the hadron mass distribution. 
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1 Introduction 



The structure of final hadronic states in semileptonic B decays 

B -» X u + I + v, (1) 

where X u is any hadronic final state coming from the fragmentation of the up quark and I + v is a lepton- 
neutrino pair, is usually studied by measuring different spectra, such as the charged-lepton energy distribution, 
the dilepton invariant mass spectrum, the hadron mass distribution, etc. The computation of these spectra is 
often rather involved because of the occurrence of different long-distance effects, both perturbative and non- 
perturbative, such as large logarithms, Fermi motion, hadronization, etc. These effects are often substantial in 
the end-point regions relevant for the experimental analysis. 

Let us summarize in physical terms the main properties of semi-inclusive B decays E] ■ We begin by 
considering the radiative decays 

B -» X s + 7 (2) 

which have a simpler long-distance structure than the semileptonic decays Q . That is because, in the radiative 
case, the tree-level process is the two-body decay 

b s + 7 (3) 

having the large final hadron energy 

2 E x = m b (l - ^ S* m b , (4) 

where in the last member we have neglected the small (compared with the beauty mass mj) strange mass m s . 
The final hadronic energy Ex fixes the hard scale Q in the decay: 

Q = 2 E x , (5) 

so that Q = m b . We are interested in the threshold region 

m x < E x , (6) 

which can be considered a kind of "perturbation" of the tree-level process due to soft-gluon effects, where mx 
is the final hadron mass. The hadron mass mx, which vanishes in lowest order, remains indeed small in higher 
orders because of the condition © , while the large tree- level hadron energy is only mildly increased by soft 
emissions. The perturbative expansion is therefore controlled by the coupling 

a(Q) = a(m b ) S 0.22 < 1, (7) 

where a — as is the QCD coupling, and it is therefore legitimate. 

In semileptonic decays, the tree-level process is instead the three-body decay 

b -> u + I + v (8) 

and the hadron energy Ex (i.e. the energy of the up quark) can become substantially smaller than half of the 
beauty mass m&/2. In other words, kinematical configurations are possible with 

E X « ^ (9) 

as well as with 4 

E x < m b . (11) 

4 It is clear that the condition 

E x -> A (10) 
must always be verified in order to deal with a hard process and to be justified in the use of perturbation theory. Since the beauty 
mass rrif, is larger than the hadronic scale A by an order of magnitude only, it is in practise not easy to satisfy all the conditions 
above. In any case, the theory is constructed by taking the limit mj — * oo. 
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Consider for example the kinematical configuration, in the b rest frame, with the electron and the neutrino 
parallel to each other, for which condition © holds true, or the configuration with the leptons back to back, 
each one with an energy sw m&/2, for which condition is instead verified. This fact is the basic additional 
complication in going from the radiative decays to the semileptonic ones: the hard scale is no more fixed by the 
heavy flavor mass but it depends on the kinematics according to eq. JSJ . 

Before going on, let us give the definition of a short-distance quantity in perturbation theory. Let us consider 
a process characterized by a hard scale 

Q > A, (12) 
where A is the QCD scale, and by the infrared scales 

^F, rn\ « Q 2 (13) 

— the soft scale and the collinear scale respectively. If infrared effects are absent in the quantity under consid- 
eration, logarithmic terms of the form 

O 2 

a"log fc -V (n = 1,2, ••• oo, k = 1,2, ••• 2n) (14) 

m x 

must not appear in its perturbative expansion. The presence of terms of the form (|14|) would indeed signal 
significant contributions from small momentum scales. That is because these terms originate from the integration 
of the infrared-enhanced pieces of the QCD matrix elements from the hard scale down to one of the infrared 
scales: 

dk 2 ± [ Q2 dk\ , Q 2 

If perturbation theory shows a significant contribution from small momentum scales to some cross section or 
decay width, we believe there is no reason to think that the same should not occur in a non-perturbative 
computation. 

On the contrary, a quantity such as a spectrum or a ratio of spectra is long-distance dominated if it contains 
infrared logarithmic terms in the perturbative expansion of the form l|14fl • These terms must be resummed to 
all orders of perturbation theory in the threshold region ijBJ . For hadronic masses as small as 

m x « VAQ, (16) 

the non-perturbative effects related to the soft interations of the b quark in the B meson (the so-called Fermi 
motion) also come into play and have to be included by means of a structure function 5 . Our criterion to 
establish whether a quantity is short-distance or it is not is rather "narrow" ; we believe that it is a fundamental 
one and we use it sistematically, i.e. we derive all its consequences. The consequences, as we are going to show, 
are in some cases not trivial. 

Factorization and resummation of threshold logarithms in the radiative case is similar to that of shape 
variables in e + e~ annihilations to hadrons. One can write: 

1 f*' dT 



, ,J-dt' s = C r [a(m b )} £ [t s ; a(m b )] + D r [t s ; a(m b )} , (17) 

1 T JO Ut s 



where T r is the inclusive radiative width, we have defined 

o 

m 



(0 < t. < 1) (18) 



m 2 



and: 



5 The first and the second integral on the l.h.s. of eq. H15I are related to the structure function and to its coefficient function 
respectively. 
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• C r (a) is a short-distance, process-dependent coefficient function, independent on the hadron variable t s 
and having an expansion in powers of a: 



C r (a) = 1 + ^Cl"'a" (19) 



n=l 

The explicit expression of the first-order correction C r has been given in 00]; 
• a) is the universal QCD form factor for heavy flavor decays and has a double expansion of the form: 



oo 2n 1 

E[«;a] = 1 + £ ^ a" log" - 



n=l k=l 



1 aC F , 2 1 7aC Fl 1 1 /aC F \ 2 , 4 1 



= 1-x -log^- + - Mog- + - M log 4 - + -.., (20) 

2 7T U 4 7T U 8 \ 7T / U 

where CV = (iV^ — l)/(2Nc) is the Casimir of the fundamental representation of the group SU(3) and 
Nq = 3 is the number of colors. In higher orders, as it is well known, £ contains at most two logarithms 
for each power of a, coming from the overlap of the soft and the collinear region in each emission; 

• D r (t s ; a) is a short-distance-dominated, process-dependent remainder function, not containing infrared 
logarithms and vanishing for t s — > as well as for a — > 0: 

oo 

O r (t,;a) = 5i]jW(t.)a n . (21) 



n=l 



The first-order correction D r (t s ) has been explicitly given in 00]. 

Factorization and resummation of threshold logarithms in the semileptonic case is conveniently made start- 
ing with distributions not integrated over Ex, i.e. not integrated over the hard scale Q. The most general 
distribution in process is a triple distribution, which has a resummed expression of the form ^[2113]: 

— / — — - — — du = C \x, w] ct(w mi,)] S kt; ot(w m-b)] + D [x, u, w; a(w rrib)] , (22) 
r J Q dxdwdu' 

where we have defined the following kinematical variables: 



» = ^(0<*<1); W =^(0< W <2); u = j^V^^^ g ( U*)' (p < u < 1} . (23 ) 

mb 1 + yjl - (2m x /Q) 2 V Q J 

In the last member we have kept the leading term in the threshold region mx *C Q only. Q, the hard scale, is 
given by eq. |J5J. T is the total semileptonic width, 



= J-o 



a Cf f 25 7r i , 



(24) 



where To = G F ml\V u b\ 2 / (192ir 3 ) is the tree-level width and m& is the pole mass. Eq. (1221 is a "kinematical" 
generalization of eq. (|17JI together with the quantities involved, which all have an expansion in powers of a: 

• C [x, w; a], a short-distance, process-dependent coefficient function, dependent on the hadron and lepton 
energies but independent on the hadron variable u. The explicit expressions of the first two orders have 
been given in PE]; 
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• £[it; a], the universal QCD form factor for heavy flavor decays, which now is evaluated for a coupling 
with a general argument Q = wnib] 

• D [%, u, w; a], a short-distance dominated, process-dependent remainder function, depending on all the 
kinematical variables, not containing infrared logarithms and vanishing for u — > as well as for a — > 0. 
The first-order correction has been explicitly given in 

The resummed formula on the r.h.s. of eq. ill'L'll has been derived with general arguments in an effective field 
theory which hold to any order in a; one basically considers the infinite-mass limit for the beauty quark 

nib — > oo, (25) 

while keeping the hadronic energy Ex and the hadronic mass mx fixed. The main result is that the hard scale 
Q enters: 

1. the argument of the infrared logarithms factorized in the QCD form factor a], 

L ee logi S log -Si; (26) 
u m x 

2. the argument of the running coupling a = a(Q), from which the form factor S[u; a(Q)] depends (as well 
as the coefficient function and the remainder function). 

An explicit check of property ^ has been obtained verifying the consistency between the resummed expression 
on the r.h.s. of eq. H22|l expanded up to 0{a) and the triple distribution computed to the same order in [S]- 
Since the dependence of the coupling on the scale is a second-order effect, point [5] cannot be verified with the 
above computation. The possibility for instance that the hard scale Q — w nib in (|22f> is fixed instead by the 
beauty mass mb, 

Q = m b (???), (27) 
cannot be explicitly ruled out by comparing with the 0(a) triple distribution, because 

a{wm b ) = a(m b ) + O (a 2 ) . (28) 

A second-order computation of the triple distribution is not available at present and therefore a direct check of 
2. is not possible. An indirect check can however be obtained as explained in the following. Any semileptonic 
spectrum can be obtained from the triple distribution in (1221 by phase-space integration. In general, different 
results are obtained for a spectrum if the coupling is evaluated at the scale or instead for example at 
the scale Q27|). It is therefore sufficient to compute a spectrum in a single variable with our resummation 
scheme and to compare with an explicit second-order calculation. As far as we known, the only available 
second-order computation is that of the 0(a 2 nj) corrections to the distribution in the light-cone momentum 
P+ To compute the fermonic corrections, one simply inserts in the gluon lines of the first-order diagrams 
a fermion bubble, avoiding therefore the computation of the complicated two-loop topologies. It is however not 
straightforward to extract information about resummation from a fermionic computation, as the latter is not 
consistent "at face value" with the exponentiation property of the form factors 6 . 

In sec.|3we resum the spectrum inp + in next-to-leading order (NLO), which involves an integration over the 
hard scale Q. We also give the coefficients of the large logarithms a n \og k m,b/p+ up to third order included. The 
only coefficient that we are not able to compute is that of the single logarithm in third order, a 3 \ognrib/p+, for 
which a three-loop computation is required. These logarithmic terms can be used in phenomcnological analyses 
which do not implement the full resummed formulas or can be compared with fixed-order computations as soon 

6 Exponentiation in QED originates from the independence of mutiple soft-photon emissions. In QCD, multiple gluon emissions 
are instad not independent because of non-zero gluon color charge, but exponentiation still holds because of the cancellation of the 
non-abelian correlation effects related to the inclusive character of the form factors. 
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as the latter become available. By comparing the fermionic contributions to the coefficients of the infrared 
logarithms with the result of the Feynman diagram computation, we obtain an explicit and non-trivial check of 
our resummation scheme. 

In sec. we resum to NLO the charged lepton energy spectrum. Also in this case there is an integration 
over the hard scale Q and many considerations made for the p+ distributions can be repeated for this case as 
well. The electron spectrum is one of the first quantities studied in semileptonic decays both theoretically and 
experimentally ^JElEj; because it allows for a determination of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix element \V u b\. To avoid the large backgroud from the decay 

B X c + I + v, (29) 

it is usually necessary to restrict the analysis to electron energies above the threshold of the previous process 
— recently somehow below — i.e. to lepton energies in the interval 

The actual hadron kinematics is obtained with the replacements m c — > mo and m\, — ► ms- That gives a rather 
tight window in the end-point region, 

AE e = -» = 330 MeV, (31) 

which is dominated by large logarithms and Fermi-motion effects. We will show however that these long- 
distance phenomena are expected to have a smaller effect on this distribution than for example on the hadron 
mass spectrum. 

In sees. |2 and 13 we also discuss the validity of the so-called Brodsky-Lepage-Mackenzie (BLM) scheme jTI] 
for the evaluation of spectra in the threshold region — see also ^21 for a discussion on the radiative case @ • In 
general, this scheme aims at an estimate of the full second-order corrections 0(a 2 ) by means of the evaluation 
of the fermionic corrections, i.e. of the 0(a 2 nf) terms only. The idea is that the fermionic corrections carry 
with them a factor /3o, so after the fermionic computation, one makes the replacement: 

3 33 
n f - --(47TA,) =n f - y , (32) 

where 47r/?o = (11/3 Ca — 4/3Tr71./) is the first coefficient of the /3- function. Ca — Nc — 3 is the Casimir 
of the adjoint representation of the color group SU(3)c and = 1/2 is the trace normalization of the 
fundamental generators. The BLM approximation is often a rather good one in estimating total cross section 
and inclusive decay widths. The total semileptonic decay width for example is known at present to full order 
a 2 : 

— = K(a) = 1 + K {1) a + (n f K {2) + SK (2 A a 2 . (33) 

The explicit value of K^> has been given in eq. JUJ and the values of and SK^ can be found in |13) . 
With the BLM ansatz, is over-estimated by « 15% with respect to the exact result (with nt = 4), i.e. 
this approximation works pretty well in this case. 

In sec. 0] we study the BLM ansatz for the radiative and semileptonic distributions in the hadron invariant 
mass and for the semileptonic distribution in the variable u defined above. The conclusion is that The BLM 
ansatz works in general rather well, typically within an accuracy of 25%. 

Finally, in sec. we draw our conclusions and we discuss the forthcoming application of our formalism to 
phenomenology. 



5 



2 p + distribution 

Recently, the spectrum has been computed in the normalized light-cone momentum p + |10j defined as 

E x - \px\ uw 

P+ = = — ~ uw (0 < p+ < 1), (34) 

where in the last member we have kept the leading term in the threshold region u <C 1 only. Even though this 
variable is not special in the present framework, let us discuss its resummation for the reasons discussed in the 
introduction. For notational simplicity let us make everywhere in this section the replacement p + — > p. The 
derivation is similar to that of the resummed hadron mass distribution made in The distribution in p is 
obtained integrating the distribution in the hadronic variables u and w: 

i dr r 2 , r 1 , i d 2 r ( uw 



r dp+ J dw i max(0 ^-i) du r dudw 5 v l + u)' 

For techinical reasons, it is actually more convenient to compute the partially-integrated distribution given by: 

1 d 2 T 



1 f p dT f 
Rp{p) = p J -j-, d P' = J dudw 



r du dw V l + u 



(36) 



Let us replace the resummed expression for the distribution in the hadronic variables u and w |2j on the r.h.s. 
of eq. l|3l)l) : 



f 2 f 1 I a ir 

Rp(p) = / dw / du Ch(w; a) a [u, a(w mb)] 9 ( p 

JO imax(O.io-l) 



l + u 



2 r l / 

I I WW 

dw / dudn{u,w; a) 9 Ip — ), (37) 

imax(0,to-l) V 1 + U , 



where: 



• Ch(u>] a) is a short-distance coefficient function having the following expansion in powers of a: 

C H (w; a) = C { h\w) + aC^iw) + a 2 cff (w) + 0{a 3 ), (38) 

with 

C^iw) = 2w 2 (3~2w); (39) 



Cg\w) = —2w 2 ('S-2w) 

7T 



Li 2 (w) + logwlog(l -w) - — - - — — logw 

8 6 — 4w 



(40) 



L'i2(w) = Y^=i w n /n 2 is the standard dilogarithm; 

• a [u; a] is the differential QCD form factor: 

a [u, a] = 4- "] = S(u) + 0(a); (41) 
du 

• dH{u,w] a) is a short-distance remainder function whose expansion starts at 0{a) and which has an 
integrable singularity for u — > 0. The first-order correction for example has a logarithmic singularity: 

d^(u,w) w logu for u ^ 0. (42) 



7 The explicit expression can be found in 



G 



We are interested in the region 

p < 1. (43) 
The integral of the remainder function on the r.h.s. of eq. I|37|) vanishes in the limit p — > 0, because: 



dw du dfj(u, w: a) 6 p — / dw du dn(u, w: a) lim dip = 0. 

v 1 V 1 + u J J ' p-o V 1 + u J 



(44) 



The limit p — > can be taken inside the integral because the integrand has, as already said, at most an integrable 
singularity for u — > 0. The first integral on the r.h.s. of eq. 1|37[) can be decomposed as: 

I dw / duCniw; a) a [u, a(w nib)] 9 [ p )+ du dwCH(w;a)a[u,a(wmb)]6 ( p 

a Jo \ 1 + uJ Jo Ji V 1 



u 

(45) 

The second integral vanishes in the limit p — ► because it extends to an infinitesimal domain in w. The 
lowest-order term, obtained replacing for a[u; a] the last member of eq. (|41|1 . identically vanishes. 

To summarize, as far as logarithmically enhanced terms and constants terms in the limit p — > are concerned, 
we can make the following approximations: 8 

1. neglect the remainder function dH(u,w; a), whose contribution will be includere later on together with 
contributions of similar size; 

2. approximate the integration domain with a unit square: 

< u, w < 1; (46) 

3. simplify the kinematical constraint according to the last member in eq. (|34|l . since large logarithms of p 
can only come from the region u< 1. 

We then obtain: 

Rp(p; a) — dwCniw; a) + dwCn{w] a)T,[p/w; a(wrrib)] + 0(p; a), (47) 
Jo Jp 

where 0(p; a) denote non-logarithmic, small terms to be included later on by matching with the fixed-order 
distribution. To the same approximation, i.e. up to terms which 0(p; a), we can further simplify the distribution 
by neglecting the first integral and integrating the second integrand down to w = 0: 

Rp(p; a) = I dwCniw; a) S [p/w; a(wm,b)] + 0(p; a), (48) 
Jo 

where a = a(mb). The neglected term is indeed: 

dwCniw; a) < 1 — S [p/w; a(wmb)] f. (49) 

Expanding the form factor in powers of a(nib), the above expression reduces to a linear combination of integrals 
of the form 

rp 

/ dwC H (w; a) a n log k (p) \og l (w) ^0 for p -» (n > 1; k, I > 0). (50) 

Jo 

We have therefore shown that the neglected term is 0(p; a). 

8 It is remarkable that the variable p keeps unitary range after the simultaneous approximations 2. and 3. 
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By inscring the first-order expressions for the coefficient function Ch{w, a) and the universal QCD form factor 
a), we obtain for the event fraction Rp(j>; a), apart from small terms 0(p; a): 



/ dw Ch (w\ a) £ [p/ 
Jo 



,r\ a{wm b )} = 1 — ( - log> + — \ogp + — ) + 0{a 2 ), (51) 



where Cf = (Nc ~ l) is the Casimir of the fundamental representation of the color group SU(3) C and 

Nc — 3 is the number of colors. Let us note that we have a different coefficient for the single logarithm with 
respect to those ones in the hadron mass distribution in the semileptonic decay Q or the radiative decay 

2.1 Minimal scheme 

We introduce a resummed form of the event fraction as: 

R P (p; a) = Cp{a) S P (p; a) + D P {p; a). (52) 
All the functions above have an expansion in powers of a: 

C P {a) = 1 + aC p 1] + a 2 C ( p ] + 0(a 3 ); (53) 
Y, P (u;a) = I + aY, { p\p) +a 2 Z P 2) (p) + 0(a 3 ); (54) 
D P (p;a) = aD ( p\p) + a 2 Df {p) + 0(a 3 ). (55) 

Let us begin considering a minimal scheme, where only logarithms in p are included in the effective form factor 
Sp. Since we will consider a different scheme later on, let us denote the quantities in the minimal scheme with 
a bar. The first-order corrections to the form factor and the coefficient function in the minimal scheme read: 

= -v(^ log2p + f l0SP 

OP = =-1.36461. (56) 

^ 7T 144 

Note that the correction to the coefficient function is rather large: for a(m,b) = 0.22 it amounts to w — 30%. 

We now expand the resummed result in powers of a and compare with the fixed order result, which is known 
to full order a )10j : 

R P {p;a) = 1 + aR { p\p) + a 2 Rf{p) + 0{a 3 ), (57) 

with 



4 1} w = - c -& 



TV 



i log 2 p + H bgp + 1 log 2 (p)p 3 (2 -p) + -L \og(p) p (4 - 23 p- Up 2 + 7 p 3 ) 



— (1 - p) 2 (463 + 246p - 240p 2 + 64 p 3 - 7 p 4 + 2p 5 ) 



(58) 



We obtain for the leading-order remainder function: 1 



7T 



log 2 (p)p 3 (2-p) + — \og{p) p {4 - 23 p- Up 2 + 7 p 3 ) 



9 To our knowledge, this is the only first-order remainder function possessing a double logarithmic term (multiplied by positive 
powers of the variable p). 
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I ^(l-p) 2 p(246-240p + 64p 2 -7p 3 + 2/) - lg p (2-p) 



(59) 



Since 

S P (1; a) = 1, (60) 

taking p = 1 in cq. (|52l) . we obtain the following relation between the coefficient function and the remainder 
function in the upper endpoint: 

C P (a) = 1 - D P (1; a). (61) 
It is a trivial matter to verify that the above relation holds true for our first-order expressions. 

A definition of the minimal scheme which does not make reference to the explicit perturbative expansion of 
Ch and S can be given as follows. Neglecting 0(p; a) terms, we have that 

l 

dw Ch{w\ ot) £ [p/w; a(wrrib)] = Cp(a)T,p (p; a), (62) 

so that the coefficient function in the minimal scheme can be defined taking p = 1 in the above equation and 
using eq. ijBUfl: 

Cp{a) = / dw Ch{w; a) £ [1/w; a(w m^)] . (63) 
Jo 

The effective form factor then reads: 

- j Q dwC H (w\a)Y>\p/w; a{wm b )} 

T,p(p; a) = . (64) 

J dw Ch(w; a) S [1/w; a{w mi,)] 

Let us comment upon the above result: 

• The universal form factor E(u; a) is evaluated for an unphysical arguement u > 1 in eq. (|64|l . but this 
does not cause any problem because the perturbative form factor is an analytic function of u which can 
be continued in principle to any (complex) value of u; 10 

• The effective form factor S p (p; a) defined in eq. 164|) has a similar form to the effective form factor for the 
t distribution £t(£; a ) considered in where t — m 2 x /rn^\ the only differences are that in the numerator 
of the definition of St, t/w 2 is replaced by p/w and in the denominator 1/w 2 is replaced by 1/w. The 
same considerations made for the ^distribution therefore can be repeated for the p-distribution. There 
are long-distance effects which cannot be extracted from the radiative decay @ , related to the integration 
over the hadron energy w in the parametric region 

p < w <C 1. (66) 

The effects of this region are however suppressed by the coefficient function: 

C H {w) w w 2 = % « 1. (67) 
mi 



10 the form factor T,(u; a) has been given to NNLO in eq.(125) of [3], where it is written as a function of 

t = f3 a(Q) log-. (65) 
u 

For u > 1 we have that r < 0, but that does not produce any singularity in S(«; a). 
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2.2 Higher orders 



Let us now discuss the higher orders. In order to remove the exponential (trivial) iterations, let us define the 
exponent of the effective form factor by means of the relation: 



with 



where we have defined: 



Wc find: 



OO n+1 



G P {p- a) = ^^<5 Pnfe a"Z£ 



71 = 1 k=l 



L p = log- > 0. 
P 



(68) 
(69) 
(70) 



Gp\2 = 

Gpn = 

Gp23 = 

Gp22 = 

Gp2\ — 



Gp34 
Gp33 

Gp32 



G\2\ 

G\\ 

G23] 

G22 

G21 
23 



12 



Ai; 



12 



z{2) 



u4 ^0o + lAi { B 1 + D l)+ A^(^-2 Z (3) 



G34', 
G33 



83 
5184 



G32 + ^A 2 (3 Q -^BiPl 



36 
Ai s 



1117 7 5 ,„, 

"20736 + 48 z(2) -12 z(3) 



48^ 2 + 



/ 23 



24 



0i 



— Ai ft (Bi - Di)+ Aipo— (— 
96 H v ' H 7T V432 



*(3) 



83 

1728 



A! 2 (£?!+#!)+ V A) 



47 25 , 

1 z 2 +Ai 

864 24 V V 



C F / 22747 
~ \ 6912 



3z(4) 



(71) 
(72) 
(73) 
(74) 



(75) 
(76) 
(77) 



(78) 



The Gij 's are the coefficients of the threshold logarithms in the exponent of the form factor G for the radiative 
decay ©• They have been explicitly given, together with all the constants Ai, Bi, Di and Pi entering the 
above equations, in z(s) = Y^=i l/ nS 1S ^he Riemann Zeta-function with z(2) = tt 2 /Q = 1.64493- • •, 

z(3) = 1.20206- •• and z(4) = tt 4 /90 = 1.08232 

Replacing explicit values for the coefficients Gpij 7^ Gij, we obtain: 



Gpn — 

Gp22 = 



C F 13 

7T 6 

Cf_ 
2tt 2 



Ca 



25 z{2) 



24 



G 



P21 



2tt 2 



V 144 4 V ' 2 



) +c ,( 



4 



1303 17 „ ,A /113 2(2) 



(79) 
(80) 

;(8i) 



llr The next-to-next-to-leading-order (NNLO) coefficients A3, B2 and D2, which are a necessary input for the determination of 
the Gjj's, have been computed in the past few years in refs. 1171 . 1181 and II 91 respectively. 
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Gp33 — 



7i~ 3 



n 2 f ( 77 11 \ / 29 

6& +C * C * 1152- g^ +C 'M 576 



/185 z(2)\ 2 / 1589 11 . 2 
C ^ n / Tone To - ) A ( -^7^ + TTa z (2) + C F 



( - r P32 — — r 



V 1296 12 J A \ 2592 24 
359 2(2) \ 2 / 64331 313 



2592 36 



\ 2 „ ^ /64331 313 , n , 11 5 , A 



(2)\ 




4 ) + 




83 


z(3) 


5184 + 


3 




5 




+ 4 



(82) 



/ 3829 2 . , 5 .A ^ 2 / 65381 7 13 11 , A . 

+ G, B/ f-H + ^( 2 ) _ f(3) V ^(^34907 _ 293 77 _ 11 

; ^ 648 144 u 24 J A V 10368 288 w 48 w 16 w 

Let us make a few remarks: 



(83) 



• As in the case of the distributions in the variables u and t studied in the leading logarithms a n L n+1 
have the same coefficients as in the hadron mass spectrum in the radiative decay l|2)l: differences only 
occur in NLO and beyond; 

• The differences between the coefficients Gpij and the coefficients G,j of the radiative spectrum are much 
smaller than for the u or the t distribution resummed in That is an indication that the long-distance 
effects in the p distribution and in the radiative decay J2J are not very different ^U] ■ 

2.3 Check of Resummation Formula 

In this section we check our results for the resummed p spectrum against an explicit computation of the fermionic 
corrections, i.e. of the 0(a 2 rif) contributions |10j . The computation with Feynman diagrams of the 0(a 2 rif) 
terms involves: 

1. neglecting double emissions off the primary charges, which would give contributions proportional to a 2 
without any n/ factor. The expansion up to 0(a 2 ) of the r.h.s. of eq. (|68|l can be made by expanding it 
in powers of Gp, 

E P = 1 + G P + ^Gp + ••• (84) 

and keeping terms up to second order in a in Gp and up to first order in Gp. There is not any multiple 
emission left in the fermonic computation, because one neglects the terms Gj, 12 , G P11 — 0(n®) in eq. 
That means that relation <|68|l is truncated to its (trivial) first order: 

Ep ~ 1 + G P [0(a 2 rif) computation] ; (85) 

2. taking into account the secondary branching of a gluon into a (real or virtual) quark-antiquark pair, 

9* -> q + q, (86) 
while neglecting the (non-abelian) gluon splitting, 

9* ^ 9 + 9, (87) 
Only a single secondary branching of abelian kind l|86|l is therefore allowed. 
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According to eq. H85|l , the best interpretation of a fermionic computation is that of relating this result directly 
to Gp and not to Ep. 

Because of|5|, a non-trivial strucure of G emerges: 

G P = G P12 aL 2 p + G P11 aL p + G P23 a 2 L 3 p + G P22 a 2 L 2 p + G P21 a 2 L p . (88) 

An 0(a 2 nf) computation therefore does not provide any check of the exponentiation property, but it gives some 
information about the structure of G. The exponentiation property of form factors has already been checked 
with many second-order and third-order computations and it is not our concern. We are only interested in the 
structure of G which is, as we have already said, sensitive to the specific integration over the hard scale. 

Our second-order coefficients Gp 2 i have a contribution proportional to n/ which is in agreement with the 



explicit computation: 

Gp 23 = + (89) 

G P22 = -g^a«/ + O(n f ); (90) 

= {h-m?) n > + °W> (91) 



where the numerical values of Gp = 4/3 and Ca = 3 have been replaced. If we assume instead that the hard 
scale Q in the coupling a = a{Q) entering the resummed formula ill'L'l) is equal to the heavy flavor mass m,b 
and not to (two-times) the hadron energy 2Ex, we obtain different fermionic contributions to Gp22 and Gp2i 
which are in disagreemnt with the fixed-order computation: 

G P t mb) = -2^»/ + 0(»J); (92) 

This provides us with an explicit check of the resummation formula for the triple-differential distribution at 
the two-loop level. Let us remark that this check is not related to the study of the BLM scheme: we just 
compared the coefficients of two independent 0(a 2 rif) computations and we never made the BLM ansatz. Our 
aim indeed was simply that of studying the effects of eq. (J5J in a partial contribution, not that of estimating 
full second-order corrections. 



2.4 The BLM scheme for the form factor 

In this section we compare our exact results for the resummed spectrum in p with those ones coming from the 
BLM ansatz. Let us make a comparison of the coefficients Gp2i in QCD in the minimal scheme defined above 
with the BLM ones Gp2i starting from the leading terms : 



Gp 23 = -t^t + tt^j = -0.185756 + 0.0112579 n f (94) 



1. a 2 L p ; the QCD coefficient 

11 + JtL 
6tt 2 9tt 2 

is in complete agreement with the BLM estimate: 

Gp23 = Gp23- (95) 

The approximation turns out to be exact because the QCD coefficient is proportional, as it is well known 
|20| . to the first coefficient of the (3- function (3 ; 
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2. a 2 L 2 ; the QCD coefficient is 



1 239 nr 

&P22 = 77 + tttt^ ~ tt4 = 0.242739 - 0.0168869 n f (96) 

54 108 7T 2 6 71^ 



to be compared to the BLM estimate: 

1 / 33\ 

Gp 22 = -7- [n f - y J = 0.278633 - 0.0168869^. (97) 

The BLM ansatz over-estimates the non-abelian contribution by « 15%, i.e. it works rather well; 
3. a 2 L p ; in QC-D we have the value: 

°™ -H + S - + (* - is?) - — - — ' « 

to be compared to the BLM estimate 

/ 1 11*5 \ / 33\ 

Gp2i = (^27 - TJ^J - yj = 1-13809 - 0.068974971/. (99) 

The BLM ansatz works with an accuracy better than 10%. Let us note also that the BLM cannot 
reproduce the terms proportional to the trascendental constant z(3) defined in the previous section. 

We conclude that the BLM scheme works rather well for the p spectrum. 
2.5 Non minimal scheme 

In this section we consider a non-minimal factorization scheme which has various advantages. The event fraction 
has a resummed form equal to the one in the minimal scheme, 

R P (p; a) = C P {a) £ P (p; a) + D P (p- a), (100) 

with different expressions for the coefficient function, the form factor and the remainder function. As in the 
minimal scheme, all these function have a perturbative expansion in a: 

C P (a) = 1 + aC p 1] + a 2 C { p ] + O (a 3 ) ; (101) 
Sp(p; a) = 1 + aS£ } (p) + a 2 Y,f{p) + O (a 3 ) ; (102) 
Dp (p; a) = a £>£> (p) + a 2 (p) + O (a 3 ) . (103) 

The coefficient function is defined, to all orders, as: 

C P (a) = [ dwC H (w;a), (104) 
Jo 

and the expression for the effective form factor is: 

Jo dwC H (w; a) £ [p/w: a(wm b )} 

(p; a) = . (105) 

J ™ C ff (w; a) 
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S [it; a] is the extended form factor defined in y], equal to the standard one £ [u; a] for argument less than one, 
u < 1, and equal to one for a larger argument, u > 1. Explicitly (cfr. eq. (07))) 12 : 



E P (p; a) 



dwCn (w; a) + dwCn (w; a) £ [p/u>; afwmi,)] 



/ dw (Tff (iy; a) 

This effective form factor is normalized, to all orders, as in the minimal scheme: 

E P (1; a) = 1. 



(106) 



(107) 



By inserting the expansion for Ch(w; a) and for S(u; a) and integrating over the hadronic energy w, we obtain 
for the first-order terms: 



C« = - — 



CV 335 



7T 144 



and 



4%) 



-0.98735 

CWl , 2 13 , 8 25 , 

- log zH log n H p 

it \ 2 6 ^ 9 18^ 



(108) 



(109) 



The coefficient function has a smaller value in this scheme than in the minimal one and the effective form factor 
contains not only a constant term for p — > but also infinitesimal terms in the same limit. Matching with the 
fixed-order distribution as we have done in the minimal case, we obtain for the remainder function: 



d p (P) = 

7T 



- ^(2-p)p 3 log 2 p- -L(4-23p-14p 2 + 7p 3 )plo^, 



(680 + 269 p - 990 p 2 + 447 p 3 - 80 p 4 + lip 5 - 2p 6 ) 



(110) 



Let us note that, with our scheme choices, 

C P (a) = C T (a) = Cu(a) (111) 

to all orders in a, where Ct{oc) and Cjj{ol) are the coefficient functions for the distributions in the variables u 
and t computed in 

The exponent of the effective form factor in this non-minimal scheme Gp contains, in the limit p — > 0, also 
constants terms (Gpij — Gpij for j > 1): 



23 5 

Gpio = - — A x + — (B X + A); (112) 

23 101 23 / 2 \ 47 7 

Gp 2 „ = - Tii ^-^^A + ^/3o(i? 1 + 3^j^^(B 1+ ^) + ^(5 1 + ^) 2 + 

+ ^ (Ba + £ 2 ) + ^ ,( 2 ) ^ (Bx + Bx) + Bx ^ - 2 2 ( 3)) + A ^ - 2 ,(3)) + 



Explicitly: 



t 2 / 2057 23 ■> ,.,,\ , f /' 901 21 5 



z(2) + — 2(3) + Ax — + - 2(3) + 3 2(4) . (113) 

V 41472 144 w 12 v 7 tt I 20736 6 v ' w ' 



G P10 = (114) 



12 Rcplacing the extended form factor £ with the usual one £ in eq. I1U5I amounts to dropping the infinitesimal terms for p — » 
in Ep. 
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Gp20 



C F 



/2243 5 \ _ / 14435 
»/ 1^7" ™*( 2 ) )+ C F — 



V5184 72 



1728 



83 , , 19 , s , N 5 , A 
— z(2) + T z(3) + 3z(4)-- 2 (3)j + 



+ C A 



24775 



193 5 ,„. 

10368 ' 288 Z(2) + 48 Z(3) 



(115) 



As in the case of the t distribution considered in 0, the coefficients functions are related in the two schemes 
by the following relations: 



Op — Op — Opio, 



,(i) 



r {2) _ M2) 



(116) 
(117) 



The first of the above equations is easily verified by inserting our first-order expressions. 



3 Electron spectrum 

In our framework, the electron spectrum is obtained by integrating the distribution in the hadron and electron 
energies w and x = 1 — x, resummed to NLO in [2], over the hadron energy: 

ldT = f 1 + i dw ]_^£_ (n8) 

r dx J s r dwdx 

In the previous case, in order to avoid generalized functions entering the differential form factor a[u; a], we have 
considered the partially integrated distribution or event fraction. That in not necessary in this case because 
the differential electron spectrum, as we are going to show, already contains the partially integrated form factor 
£[m; a]. Inserting the resummed expression for the double distribution, we obtain: 13 



1 dT 

r dx 



I dw CL.(x,w;a)T, —;a(wmb) + / dw Cxwi(x; a)\l — Cxw2{x; a)S [w — 1; a(rrib)] \ 



+ / dw d < (x,w; a) + / dw d > (x,w; a). (119) 
We are interested in the region 

X < 1, (120) 

which selects hadron final states with a small invariant mass and produces large logarithms in the perturbative 
expansion. The integral of the first remainder function d<(x, w; a) vanishes in the limit x — > 0: 

i ,i 
dwd < (x,w;a) ~ / dwd < (x,w;a) — » for x — * 0, (121) 



Jo 
because: 

d<(x,w;a)^0 for x -► (0 < w < 1). (122) 

The integral of the second remainder function d > (x 1 w; a) also vanishes in the same limit: 

1+2 

dw d > (x,w; a) — > for x — > 0, (123) 



13 The exploit expressions of the coefficients functions and the remainder functions can be found in 
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because the integral extends to an infinitesimal region and the integrand vanishes for w — > 1 + : 

d>(x,w; a) — > for w — > 1 + . 



(124) 



Finally, the integral involving E[u> — l;a] also vanishes because, after expanding the form factor in powers of 
a, it produces terms of the form: 



J dw a n \og k {w - 1) for x -► 0, 



(125) 



because the logarithm function to any positive power k > 0, log fe (it), has an integrable singularity in u — 0. We 
may write therefore: 



i c/r 

— — = / dw Cl{x,w; a) E — ; a(wmb) + 0(x;a), 



r dx 



. w 



(126) 



where by 0(x; a) we denote terms which vanish for x — > as well as for a — > (the neglected terms also vanish 
for a — ► 0). Let us conclude that, as long as logarithmically enhanced terms for x — > and constants are 
concerned, the integration can be made in the leading region 



< w < 1, 

which coincides with the phase-space domain at the tree level. 



(127) 



3.1 Minimal scheme 



For clarity's sake let us consider at first a minimal factorization scheme. Since we are in region l|120[l . we can 
take the limit x — > in the r.h.s. of eq. I|126|) anytime we do not encounter singularities: 



1 dT 1 f \x 

— — = - I dw Cl(0, w; a) S — ; a{wm b ) + 0(x), 
21 dx 2 J lw 



(128) 



where we have divided by a factor two in order to simplify the forthcoming formulas. Let us note that, since 
we have taken the limit x — ► in the coefficient function Cl{x, w; a) which also has an 0(a ) contribution and 
we have integated down to w — 0, we have modified the integral above by terms 0(x) not multiplied by a: we 
will take care of them in a moment. Replacing the first-order expressions for the coefficient function [3], 



C L (0,w; a) = 12w(l - w){ 1 + 



T ■ , \ i f-, \ ^ wlogw 35 

Li 2 (w + log w log 1 - w - -logw - — - — 

2 2(1 — w; ) 8 



+ O (a 2 ) 



and the QCD form factor given in eq. (|20() . we obtain: 



Cl(0, w; a) E — ;a(wmb) 
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as Cf { 1 , i 31 , 15\ 
= 1 - I 2 lo § x + - logx + T ) + 0(a 2 ). 



(129) 



(130) 



In order to avoid a remainder function O(a ), let us introduce an over-all factor equal to the lowest-order 
spectrum, so that the resummed form for the electron spectrum reads: 



1 dT 
2TdE 



(1 — x) 2 (1 + 2x) Cx(a)T,x(x; a) + dx(x; a) 

We require that the remainder function vanishes for x —> 0: 

lim dx(x; a) = 0. 

X— 

The coefficient function, 



C x (a) = 1 + aC ( x } + a 2 C (2 \ +0{a 3 ) 



(131) 

(132) 
(133) 



16 



has the leading correction: 



C 



(i) 
x 



15 C F 

47T 



5 

7T 



1.59155. 



(134) 



Let us note that this correction is very large as it amounts to ea —35% for a(m&) = 0.22. The effective electron 
form factor, 

(135) 



E x (5; a) = 1 + aE^{x) + a 2 Y%\x) + 0(a a ), 



has a first-order expression: 



,-u .' Cf ,2 - 31 C F 

S X W = - 7T- lo S 1 7^ ^S^- 

Z 7T 1Z 7T 



(136) 



Let us now compute the remainder function, which is necessary to describe also the region 5 ~ O(l) and to 
have a uniform approximation in the whole 5 domain. The remainder function has, as usual, an expansion of 
the form: 

^ (137) 



d x (x; a) = ad^ } (x) + a 2 S 2) (5) + 0{a A ). 



We match with the fixed-order spectrum, 



^ J = (1-5) 2 (1 + 25) [l + a s H^{x) + a 2 H^(x) + 0(a 3 ) 



whose first order correction is well known E] : 



(138) 



7T 



log x 



41 



41 



36(l-5) 2 
23 



36(1-5) 18(1 + 2 5) 



13 



55 



54(1-5) 54(1 + 2 5) 
• log 5 log(l - 5) + Li 2 (5) 



log a; + 



(139) 



The leading-order remainder function reads: 



— < log 5 log(l - 5) + Li 2 (5) + — 

7T 1Z 



41 



41 



13 



36(1-5) 
55 



23 t 

12 ~ 36(l-5) 2 + 54(1 - 5) ~ 54 (1 + 25) 



18(1 + 25) 
log5 }. 



+ 



(140) 



It vanishes, as requested, for x — > and, despite appearances, it has not any 1/(1 — x) singularity. 
The coefficient function can be defined to all orders as: 



1 f L 

C x (a) = - I dwC L (0,w; a)S 

2 Jo 



— ; a(w nib) 
w 



and the effective form factor as: 



T, x (x; a) 



The latter is normalized to all orders as: 



J dw Cl(0, w; a) S [x/w; a(w nib)] 
Jq dw Cl(0, w; a) £ [ 1/iy; a(w mj)] 



£*(!; a) = 1. 



(141) 



(142) 



(143) 



The effective form factor Ijl42|l factorizes the long-distance effects related to the threshold region, a part of 
which cannot be extracted from the radiative decay J2J. An integration over the hard scale Q — wnib is indeed 
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involved in Ex while, in the radiative decay, kinematics fixes w = 1 and there is no such integration. The 
relevant integration region is, parametrically: 

x < w < 1, (144) 

where the large logarithms | \og(x/w)\ S> 1 are multiplied by a coupling a(wrrib) 3> a(m&). As in the case of 
the p+ distribution, the effects of this region are however suppressed by the coefficient function: 

C L (w) « w = — < 1. (145) 



3.2 Higher orders 

Let us write as usual 

Ex = e Gx , 

where the exponent of the form factor has the perturbative expansion: 

OO 71+ 1 



G x (x; a) = ^Y,G Xnk a n Lg, 



with 



n=l k=l 



1 



L s = log - > 0. 

X 



(146) 



(147) 



(148) 



By inserting the truncated expansions for the coefficient function Cl and the form factor E in eq. 1142|) . 
expanding the product, integrating term by term and taking the logarithm, we obtain: 



Gx\2 
Gxn 

Gx23 
Gx22 

Gx21 

Gx34 
Gx33 



G 



X:',2 



G 
G 



12; 



11 



■An 



= G 
= G 



23; 



22 



T2 A * 



T2 Al(3o > 



= G 2 i + -A 2 + A\ 


r /3 i?i + ^i — 



25 



6 

= G 



6^-54 

'^-2,(3) 
216 1 ' 



34; 



^A x p + -A X (B X + D x ) 



= G 



33 



■ —A? + —A? p + — A x 
648 1 72 iro 18 Mo ' 



G32 + \a 2 (3 q - I p Q 2 B x + H A x A 



13 
72 



36 

A X P {B X -D X )+A X ^ (H--(3))+V^ (^-3*(4) 



Ai (— /3 2 + — A 
1 I 3fi MU 12 



36 

z(3) +Ai — — - 



(149) 
(150) 
(151) 
(152) 



(153) 
(154) 
(155) 



(156) 



Explicitly, one has for the coefficients Gxij 7^ Gif 
31C F 



Gxn 



12 7T ' 



(157) 
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GX22 
Gx21 

Gx33 



Cf 

IT 2 
it 3 



C F 



23 
144 



nj + Cf 



13 
72 



z(2) 



73 
72 



z(2) 



1296 
2573 



143 11 



C A 
11 



205 z(2) 
288 + ~ 
3 



+ — z(2)--z(3) +C A 



(> 



23 



11 



■z{2) 



2(3) 



(158) 
; (159) 



"/ + Cf ^l86i-T^J +C ' n n2i6 



s(2) 



, 65 z(2) 

CAn/ (648 -l2- 



Gx32 — — r 



5184 
229 



11 

24 

z{2) 



z(2) 



/j35_ £(3) 



^ U4S 



1296 36 ' "'' 
401 121 



/2807 1183 



V1296 



11 5 ,, 

288 *(2)- T *(3) + i *(4) 



1296 144 
10115 29 



5184 72 



^(2) 



2(2) 



12 

2(3) 
24 



2(3) 



Cf \288 

elf 6217 



13 
f 36 

29 



2(2) 



31 
12 

77 



2(3) 



11 



2(4) 



V 1296 24 



z (2) + -z(3) 



11 

16 



2(4) 



(160) 



(161) 



The coefficient of the single logarithm at 0(a), Gxu, is the same as in the distribution in the hadron mass 
squared t, while the higher-order coefficients are different. Therefore there is not a simple relation between 
these two spectra beyond leading order from two loops on, as was instead guessed in |22| . 



3.3 Study of the BLM scheme 

In this section we compare the exact second-order coefficients in the exponent of the electron form factor Gx 
with those ones obtained by the BLM ansatz: 

2053 1 23 , 
G X 22 = 648^ + 54~108^ n ' = 0.339525 - 0.0215777^ (162) 

23 / 33 \ 

G X 22 = -—^2 U/ - yj = 0.356031 - 0.0215777n/ (163) 

G X 2i = - (j^g - 5f) ( n f ~ y) = L64891 ~ 0.0999342^ (165) 
As in the previous case, this approximation works within O(10%), i.e. rather well. 



3.4 Non-minimal scheme 



As discussed in detail in j2j, for phenomenological applications it is convenient to dehne a non- minimal scheme 
which involves the integration of the universal form factor T,(u;a) only in the physical region < u < 1. Let 
us introduce also for this scheme an over-all factor equal to the spectrum in lowest order: 



I^=(l 
2T dx K 



x) 2 (1 + 2x) Cx{ot) Yix(x; a) + dx(x; a) 



We choose to dehne the coefficient function as: 

C x (a) e 



dw Cl(0, w; a). 



(166) 



(167) 



19 



One immediately obtains for the first-order correction: 



Cj? = -— ^ = -0.748618. (168) 

7T 72 



The correction to the coefficient function is — as in all the cases we have considered — negative. For a(mb) 
0.22 the correction is « —16.5%, i.e. it is basically half of that in the minimal scheme. 

The effective electron form factor can be defined as: 
Jq dw Cl{0,w; Qi)S[x/io; a(wmb)] 



T, x (x; a) — , 

J Q dwC L (0,w; a) 

f* dwC L (0,w; a) + / $ dwC L (0,w; a)T,[x/w; a(wm b )] 



dwC L (0,w; a) 

where a = a[rnb). In order to simplify the definition as much as possible, we have integrated down to w = 0, 
replacing the form factor a] with its extension £[u; a]. The normalization is the same as in the minimal 
scheme: 

Sjc(1; a) = 1. (170) 
The perturbative expansion of the form factor is: 

E x (x; a) = 1 + a£W(i) + a 2 £^(x) + 0(a 3 ). (171) 

By replacing the explicit expressions for the coefficient function Cl and the form factor S, expanding in a and 
performing the integration, we obtain: 

,(i)/_x CW l ra ,31 . 143 27 _ 2 25 _ 3 



7T 



E\i ; (x) = — [ - -Li + —L s + — x 2 x A . (172) 

x K ' - [ 2 * 12 72 8 18 ' 



In general, the form factor contains constant terms for x — > as well as vanishing terms in the same limit. The 
former modify the coefficient function while the latter modify the remainder function with respect to the values 
in the minimal scheme. 

The remainder function reads in this scheme: 



~m, J 29 27 2 25 3 41 23 , 

d * {x) = — 12 + 18 X - 36(1^) " 18 il + 2-x) +^xlo g (l-x)+U 2 (x) 



23 41 13 55 



12 36 (1 - x) 2 54 (1 - as) 54 (1 + 2 x) 



\ogx\. (173) 



This remainder function also vanishes for x — > and differs from the one in the minimal scheme only for the 
subtraction of the infinitesimal terms in Ijl72|l • 

The coefficients of the infrared logarithms Gxij in the exponent of the form factor Gx are the same as those 
in the minimal scheme computed in the previous section: 

G Xij = G Xij for j > 1. (174) 

In the limit x, — > (i.e. neglecting O(x) terms), Gx also contains the constant terms: 

19 5 

Gxw = ~-A 1 + - (Bi +£>!); (175) 
36 6 
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Gx2 = ~^A 2 - ^UiA, + H( Bl+Z3l ) 2 + ™p (3fl 1 + 2D 1 ) + | (B 2 + Z? 2 ) 
do 72 72 do o 
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Explicitly: 
Gxw 

Gx20 
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(177) 



(178) 



The relations between the coefficients functions in the two schemes have the same form as those for the p_| 
spectrum given in the previous section: 
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(179) 
(180) 



The first of the above equations can be directly verified by inserting the first-order quantities. 

Long-distance phenomena (large logarithms, Fermi motion, hadronization, etc.) are expected to have a 
larger effect in the distributions in the variables q = t, u or p than in the electron energy spectrum . That is 
because the former variables have a tree-level distribution consisting of a peak in zero: 



1 dT 
f ~dq 



5(q) + Q{a) 



while the electron spectrum has a broader distribution with a maximum in x = 0: 



1 dT 
2T ~dE 



(1 - xf (1 + 2x) + 0(a). 



(q = u, t, p), (181) 



(182) 



Long-distance phenomena always have a smearing effect, which is more prounanced in a starting distribution 
of the form l|181|) than of the form l|182|) . Technically, that is reflected by the fact that long-distance effects 
in the electron spectrum are tempered by the presence of the partially-integrated form factor £ instead of the 
differential one a. We therefore expect the electron spectrum to be a less sensitive quantity to threshold effects. 



4 The BLM scheme for other spectra 

The accuracy of the BLM ansatz can also be studied by comparing its predictions with the exact results 
obtained for the following distributions treated in [2] and in 0: 

14 The u spectrum have been studied in pj. 
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1. hadron mass distribution in the radiative decay J2J) (or, equivalcntly, scmilcptonic distributions not in- 
tegrated over the hadron energy). Our results for the QCD coefficients and the BLM estimates read 
respectively: 

G 22 = JL + 1 - n f = 0.152206 - 0.0121961 n /; (183) 

13 / 33\ 

G 22 = irif - y J = 0.201235 - 0.0121961 nf, (184) 

z(3) 917 35 / 85 1 
9^ + 72 7T 2 ~ 54 ~ V 108 71-2 _ 27 

= " (lbfb ~ ^7) ("/ - f) = 0-704657 - 0.0427065 n,. (186) 

The estimate of the non-abelian contribution is accurate with 0(25%); 

2. distribution in the variable u defined in the introduction (basically the hadron mass normalized to the 
hadron energy). We obtain: 



G 21 = - -12 + - - _ ) nf = 0.628757 - 0.0427065 n f ; (185) 



5 1 1 

24^2 + 54 ~ 36^2 

1 / 33 s 

36"^ [ nf ~ T 

z(3) 217 35 / 10 1 
9^2 + 12 7T 2 ~ 54 ~ V 9^ ~ 27 



Gc/ 22 = — 77-; — n + 77 - ^-o"/ = -0.00259006 - 0.00281448 n/; (187) 
G^ 22 = -tjtt^ (n/ - ^ J = 0.0464389 - 0.00281448 n /; (188) 



G[/21 = -_y + __ _ __]„ / = 1.17054 - 0.0755421 nf (189) 



Gc/21 - - (J^ - ^) (n f - f ) = 1-24644 - 0.0755421 nf (190) 

The agreement is good for the G\j%\ coefficient, within 10%, but not good for the G[/ 22 coefficient, which 
is however very small because of some accidental cancellation; 

3. hadron mass spectrum in the semileptonic case. We obtain: 



G T22 = — -5 + 77 - ^-"2 n/ = 0.514336 - 0.0309593 n/ (191) 
216 7r 2 54 36r 

11 / 33\ 

G T22 = -ttt^ [nf - — = 0.510828 - 0.0309593n /; (192) 



G T21 = - -^1 + ^ - - - (^-3 - - ) „, = 0.385075 - 0.0408302 B/; , I m , 
G ^i = - ( 77-^ - = (rif - — = 0.673698 - 0.0408302 n,. (194) 



The estimate of Gt 22 is very accurate, while for Gt 2 i it is not. 

The over-all picture is that the BLM ansatz offers an estimate of the resummation coefficients with an acceptable 
relative error most of the times. There are "exceptional cases" , when the coefficients are very small because of 
some accidental cancellation, which cannot be accounted for by the general BLM ansatz. 



5 Conclusions 

In this paper we have resummed to next-to-leading order the distributions in the light-cone momentum p + = 
Ex — \Px\ and in the charged lepton energy E e in the semileptonic decays l(T|l. where Ex and px are the total 
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energy and three- momentum of the final hadron state X u . These spectra have a different threshold structure 
with respect to the hadron energy distribution in Q or the photon spectrum in J2J), because they involve 
integration over Ex, i.e. over the hard scale Q = 2Ex of the hadronic subprocess in We have also 
presented a qualitative discussion of specific aspects of the long-distance effects in the electron energy spectrum: 
they should be suppressed with respect to the distribution in the hadron mass mx or in p + because of a smearing 
effect of the Born three-body kinematics. 

We have explicitly checked our resummation formula for the p+ form factor by expanding it to 0(a 2 ), picking 
up the nf contributions and comparing them with the 0{a 2 rif) corrections computed with Feynman diagrams. 
This is a highly non-trivial check at the two-loop level of the resummation formula for the triple differential 
distribution and makes us more confident of our formalism. 

We have also discussed the validity of the so-called Brodsky-Lepage-Mackanzie (BLM) scheme for the 
estimate of the second-order corrections to the form factors for various spectra. The conclusions are that this 
scheme offers in general a rather good approximation: it exactly reproduces the coefficients of the a 2 log 3 terms 
and gives good numerical estimates of the a 2 log 2 and a 2 log coefficients. In cases in which the QCD coefficients 
are anomaously small because of some accidental cancellation, the BLM scheme is not accurate, as in these 
cases the n /-dependent and the n /-independent contributions presumably cancel to a different degree. It would 
be interesting to extend the BLM analysis to the third-order coefficients which have been computed by us for 
all the distributions treated. Let us remark that the exponentiation property does not follow from the BLM 
ansata and it has to be imposed externally. 

By expanding the resummed formulas up to 0(a 3 ), we have computed the coefficients Gpij and Gxij 
of the infrared logarithms log nib /p+ and logl/cc in the exponent of the effective form factors Gp and Gx 
respectively. These coefficients constitute a true prediciton of our resummation schme and can be checked 
with second and third order computations as soon as they become available. These coefficients can also be 
used in phcnomcnological studies which do not implement the full resummed formulas. For both distributions 
we have defined (standard) minimal factorization schemes and also non-minimal schemes which can be easily 
implemented in phcnomcnological analyses. As in our previous work [3], we find that the perturbative expansion 
seems to have better convergence properties in the modified schemes than in the minimal one. 

Analytic expressions for double distributions or for single distributions with some kinematical cut can be 
obtained in similar way as we have made. This work completes however our study of the general factorization- 
resummation properties of semileptonic decay spectra ^|21yJ: for each distribution considered we have com- 
puted the coefficient function and the remainder function in NLO and the (effective) form factor in NNLO. The 
next step is a phenomenological analysis of the experimental data on semileptonic and radiative decays by using 
a model for the QCD form factor such as that formulated in |2.'ij . combined eventually with non-perturbative 
components not calculable in perturbatve QCD. 

Note added 



A few days after the first version of this paper was put on the archive, another work on resummation in semi- 
leptonic decays has also appeared |24|. In this paper, expressions of the second-order coefficients of the p + form 
factor are derived, which are in complete agreement with ours. 
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